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The  Principle  of  Independence  of  the  Cavity  Sections  Expansion 
(Logvinovich’s  Principle)  as  the  Basis  for  Investigation  on  Cavitation  Flows 


A.  D.  Vasin 

State  SRC  TsAGI,  Moscow,  Russia 
17,  Radio  str.,  Moscow,  Russia 


Summary 

For  the  investigation  of  complex  non-steady  cavitation  flows  in  three  dimensions  it  is  extremely 
important  to  find  a simple  model  which  on  one  hand  would  correspond  to  the  basic  laws  of  a fluid  stream 
(the  laws  of  conservation  of  energy  and  momentum)  and  on  the  other  hand  would  be  simple  enough  for 
the  numerical  calculation.  This  model  was  presented  by  G.  V.  Logvinovich,  and  was  called  “The 
Principle  of  Independence  of  the  Cavity  Sections  Expansion”.  In  the  present  lecture  the  deduction  of  the 
principle  of  independence  from  the  equation  of  conservation  of  energy  is  adduced.  It  is  shown  that  for 
the  stationary  cavities  the  results  obtained  on  the  basis  of  the  principle  of  independence  agree  with  the 
results  obtained  on  the  basis  of  the  slender  body  theory  and  the  numerical  methods  as  well  as  with  the 
experimental  data.  It  is  shown  that  for  the  special  case  - the  cavitation  number  equals  zero  - the  cavity 
shape  determined  on  the  basis  of  the  principle  of  independence  agrees  with  the  well-known  Levinson- 
Gurevich  asymptotic.  The  application  of  the  principle  of  independence  for  the  calculation  of  some  types 
of  the  non-stationary  cavities  and  cavities  with  variable  external  pressure  (the  vertical  cavities  in  the 
gravitational  field)  is  considered. 


Introduction 

In  the  case  of  high-speed  motion  in  water,  cavities  filled  with  gas  or  vapor  are  formed  behind  the  body- 
cavitator.  The  mathematical  problem  on  the  determination  of  the  cavity  shape  is  stated  as  the  inverse 
problem  on  the  flow  around  a body  - the  free  boundary  shape  is  found  from  the  constant  pressure 
condition.  In  case  of  the  non-stationary  cavities,  the  cavities  with  variable  external  pressure  along  their 
length  or  the  asymmetrical  cavities,  the  problem  on  the  determination  of  the  cavity  shape  becomes 
extremely  complex.  For  the  practical  problems  it  is  very  important  to  have  a simple  method  of 
calculation  of  non-stationary  cavitation  flows  that  corresponds  with  the  basic  laws  of  a fluid  stream. 
Such  a method  was  suggested  by  Logvinovich  and  was  called  “The  Principle  of  Independence  of  the 
Cavity  Sections  Expansion”. 

The  principle  of  independence  is  stated  as  following.  Each  cross  section  of  a cavity  expands  relatively  to 
the  trajectory  of  the  center  of  a cavitator  which  happens  almost  independently  from  the  following  or  the 
previous  body  motion.  The  expansion  occurs  according  to  the  definite  law  which  is  dependent  only  upon 
the  difference  between  the  pressure  at  infinity  and  the  pressure  within  the  cavity,  the  speed,  the  size,  and 
the  drag  of  a body  at  the  moment  when  a body  passes  the  plane  of  the  considered  section  [I].  The  law  of 
expansion  of  the  cavity  fixed  section  can  be  determined  using  the  principle  of  conservation  of  energy  in 
the  wake.  The  thin  cavity  in  ideal  fluid  can  be  considered  as  the  wake.  Each  unit  of  the  wake  length 
keeps  the  energy  which  was  expended  by  the  cavitator  at  the  moment  when  it  passes  this  unit  of  the 
length  of  the  trajectory  [2,3].  This  theory  applied  to  the  cavity  with  the  constant  difference  of  pressures 
gives  results  which  are  very  similar  to  the  results  observed  in  the  experiments.  Therefore,  we  can  suggest 
that  this  theory  is  also  applicable  for  the  case  of  variable  difference  of  pressures. 


Paper  presented  at  the  RTOAVT Lecture  Series  on  “Supercavitating  Flows  ”,  held  at  the  von  Karman 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 
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1.  The  Deduction  of  the  Principle  of  Independence  from  the  Equation  of 
Conservation  of  Energy 

We  designate  the  coordinate  along  the  arc  of  the  trajectory  of  the  cavitator  center  as  h (Fig.  1).  This 
coordinate  is  fixed  relatively  to  the  motionless  fluid.  At  the  point  h let  us  draw  the  plane  E that  is 
perpendicular  to  the  trajectory.  We  will  observe  the  development  of  the  cavity  cross  section  arisen  on 
this  plane  at  the  moment  t= 0. 


Fig.  1 The  scheme  of  development  of  the  non-stationary  cavity 

Let  us  introduce  the  following  designations:  the  radius  of  the  cavity  section  R=R(h,t),  the  area  of  the 
cavity  section  S(h,t)=nR2,  the  pressure  within  the  cavity  Pk=Pk(h,t),  the  pressure  at  infinity  P„(h,t).  We 
can  consider  the  pressure  P„(h,t)  to  be  the  pressure  at  the  point  “0”  of  the  intersection  of  the  trajectory 
with  the  plane  E when  the  cavity  is  absent  or  to  be  the  hydrostatic  pressure  on  the  depth  of  the  point  ‘O’. 

The  cavitator  passing  the  path  Ah  performs  certain  work  and  gives  the  energy  WAh  to  the  fluid  (where  W 
is  the  drag  force  acting  on  the  cavitator  at  the  moment  t= 0 ).  Proceeding  from  the  physical  scheme  of  the 
flow  we  can  accept  that  approximately  the  energy  WAh  is  conserved  in  the  considered  section  at  the  same 
segment  Ah  in  the  form  of  the  kinetic  energy  TAh  ii  potential  energy  EAh  up  to  the  moment  of  this 
section  closure.  Therefore,  we  can  write  an  approximate  equation  for  each  cavity  cross  section 


T(h,t)+E(h,t)=W(h,  0) 


(1.1) 


Equation  (1.1)  can  be  applied  to  the  whole  cavity  length  and  signifies  that  the  energy  that  was  imparted 
by  the  cavitator  at  the  point  of  the  trajectory  h is  conserved  on  the  cavity  at  this  point.  Proceeding  from 
the  first  Green’s  formula  [4]  we  can  express  the  kinetic  energy  that  fits  the  unit  of  the  cavity  length  by  the 
following  equation 


T = 


— pep  2 kR 
2 


d(p 

dn 


(1.2) 
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where  (p  is  the  potential  of  the  absolute  velocity  at  the  boundary,  R and  R~d(p/dn  are  the  radius  and 
the  radial  velocity  of  cavity  boundary  at  the  point  of  the  trajectory  h,  p is  the  fluid  density. 

The  potential  energy  of  a cavity  section  is  determined  as  follows 


E = J AP  Qi,t)2nRRdt  (1.3) 

o 


where  AP(h,t)  = P„{h,t)—  Pk{h,t) . After  substitution  of  (1.2)  and  (1.3)  in  (1.1)  we  obtain  an 
approximate  equation  of  energy 


2 


p(pS  + J A P(h,t)Sdt  = W(h, 0) 

o 


(1.4) 


where  S = 2nRR  is  the  time  derivative  of  the  area  of  a cavity  cross  section. 

The  kinetic  energy  is  determined  along  the  stream  tubes  that  lean  on  the  expanding  cavity  hole  in  the 
motionless  plane.  The  trajectory  of  a cavitator  intersects  this  plane  along  normal  at  the  point  h.  The 
dynamic  boundary  condition  is  written  in  a form  of  the  general  Bernoulli  equation  applied  to  the  points  of 
space  that  coincide  with  the  cavity  boundary  at  the  considered  moment  (later  on  we  will  be  writing  AP 
instead  of  AP(h,t)) 


d(p  v2  _ AP 
dt  2 p 


(1.5) 


where  v is  the  absolute  velocity  of  the  fluid  particles  on  the  cavity  boundary.  In  the  middle  section  of  the 
cavity  V ~ R is  the  small  quantity  and  we  neglect  the  value  of  v2  compared  to  the  value  of  A P/p.  Then 
equation  (1 .5)  is  written  in  the  form 


d(p  _ AP 
dt  p 


(1.6) 


Let  us  differentiate  equation  ( 1 .4)  with  respect  to  t and  using  expression  ( 1 .6)  we  obtain 


(1.7) 
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From  equation  ( i.6)  the  value  of  the  potential  on  the  cavity  section  boundary  is  determined  as  follows 


(P 


(1.8) 


where  (p„  is  the  value  of  the  potential  at  the  edge  of  a cavitator  at  the  moment  of  generation  of  the  cavity 
section  ( t=0 ).  After  substitution  of  (1.8)  in  (1.7)  we  obtain  the  following  differential  equation 


P<P, 


(1.9) 


t 

Taking  into  account  that  r/J  APdt  - AP dt , as  a result  of  dividing  variables  and  integrating  we  obtain 

O 


S = A 


P<Pn 


o 


(1.10) 


where  A is  the  constant  connected  with  the  initial  velocity  of  expansion  of  the  cavity  section  by 
dependence  A = Sa  / p(pn;  SQ  is  the  initial  velocity  of  expansion.  As  a result  equation  (1.10)  is 
represented  in  the  following  form 


S = S, 


, 1 AP, 

1 H dt 

<Pn  { P 


(1.11) 


The  value  of  the  potential  at  the  edge  of  a cavitator  we  can  write  as  following:  (pn  = — — aRnV( 0), 

where  a is  some  constant,  Rn  is  the  radius  of  a cavitator,  V (0)  is  the  velocity  of  a cavitator  at  the  moment 
t=0.  Since  for  the  stationary  cavity  the  flow  is  symmetrical  the  value  of  the  potential  on  the  mid-section 
equals  zero,  consequently  for  the  stationary  cavity  equation  (1.8)  can  be  written  in  the  following  form 


<Pn 


f AP  / 

- 1 — —dt 

t P 


~ — ti=~aR„V  (0) 
P 2 


(1.12) 
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where  4 is  the  time  of  expansion  of  a cavity  section  to  the  maximum  size.  In  case  of  the  stationary  cavity 
we  can  obtain  the  expression  of  the  constant  a from  equation  (1.12) 


oL, 

2 R„ 


(1.13) 


where  G = - 


2AP 

pv 2 


is  the  cavitation  number,  V is  the  velocity  of  stationary  motion  of  a cavitator. 


Lk  = 2 Vtk  is  the  cavity  length.  Let  us  determine  the  initial  velocity  of  expansion  Sa  from  the  equation 

of  energy  (1.4)  at  the  moment  t= 0.  Taking  into  account  that  W = Cx7lR~  ^ ^ \ where  Cx  is  the  drag 

coefficient  of  a cavitator  we  obtain 


S 2KCxRnV(Q)k  Sa  4nCx 


<Pn 


(1.14) 


As  a result  we  go  on  from  equations  (1.11)  and  (1.14)  to  the  well-known  mathematical  entry  expressing 
the  principle  of  independence  of  the  cavity  sections  expansion  in  the  form  of 


S = 


kA P 
P 


(1.15) 


Equation  (1.15)  with  the  coefficient  (1.14)  (the  constant  a is  slightly  dependent  on  the  cavitation  number 
and  its  value  is  selected  from  the  range  1.5+2)  is  widely  applied  for  the  investigation  of  non-stationary 
cavities  and  cavities  with  variable  external  pressure  along  their  length.  The  velocity  of  expansion  of  a 
cavity  section  is  determined  from  equation  (1.11) 


, fAP  , 
k — dt 

l P 


(1.16) 


Integrating  equation  ( 1 . 1 6)  we  obtain  the  dependence  of  the  area  of  a cavity  section  on  time 


S = S0+S0t- 


, f f AP  , , 

k j j dtdt 

0 0 P 


(1.17) 
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where  So  is  the  initial  area  of  a cavity  cross  section.  In  particular  for  the  stationary  cavity  the  value 


AP 


is  constant  and  from  (1.17),  (1.14)  and  (1.13)  we  obtain  the  well-known  result 


sk-s0 


(1.18) 


where  Sk  is  the  maximum  quantity  of  the  area  of  a cavity  cross  section  (the  area  of  the  mid-section). 

If  as  the  quantity  So  we  choose  the  quantity  of  the  area  of  the  cross  section  on  which  the  incline  of  a 
cavity  boundary  is  sufficiently  small  (slightly  deviating  from  a disk-cavitator),  then  the  calculations  from 
formula  (1.18)  agree  very  well  with  the  experimental  data  as  shown  in  Fig.  2 (in  this  figure  the  points 
show  the  experimental  data  for  the  different  cavitation  numbers). 

5-S0 


Fig.  2 The  comparison  of  the  stationary  cavity  profile  obtained  on  the  basis  of  the  principle  of 
independence  with  the  experimental  data 

Usually  the  profile  of  the  stationary  cavity  past  a disk  is  constructed  as  follows.  In  the  vicinity  of  the  disk 
the  cavity  shape  is  not  dependent  on  the  cavitation  number  and  can  be  expressed  in  the  form  [1] 


R(x)  = R„ 


' 3x^ 

1 + — 

v 


(1.19) 


where  x is  the  coordinate  measured  along  the  axis  of  symmetry  from  the  center  of  a disk.  The  cavity 
shape  determined  from  equation  (1.18)  joins  with  the  shape  expressed  by  the  dependence  R(x)  (1.19). 
We  can  consider  the  point  with  the  coordinates:  x=2R„,  R(x)  = \ .92  R„  to  be  the  point  of  contact. 
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2.  The  comparison  of  the  results  obtained  for  stationary  cavitation  flows  on  the 
basis  of  the  principle  of  independence  with  the  results  obtained  on  the  basis  of  the 
slender  body  theory  and  the  numerical  methods 

In  this  chapter  we  will  demonstrate  that  the  shape  of  the  stationary  cavity  represented  by 
expression  (1.18)  can  be  obtained  from  the  slender  body  theory.  As  an  example  we  use  the  cavitation 
flow  past  a thin  cone.  Let  us  consider  a more  general  case  - the  supercavity  past  a thin  cone  in  subsonic 
compressible  fluid  flow.  We  apply  the  Riabouchinsky  scheme  (Fig.  3,  the  cavity  is  closed  by  a body  with 
the  same  dimensions  as  the  caivtator). 


Fig.  3 The  Riabouchinsky  scheme  for  the  cavitation  flow  past  a thin  cone 


The  origin  of  the  orthogonal  system  x,  r is  placed  at  the  apex  of  the  cone,  the  geometric  dimensions  are 
scaled  by  the  radius  of  the  cone  base  (the  radius  of  the  cone  base  R„  equals  unity,  / is  the  altitude  of  the 
cone,  a is  the  apex  half-angle,  Lk  is  the  cavity  length,  l+Lk  and  21+Lk  are  the  coordinates  of  the  base  and 
the  apex  of  the  closing  cone  respectively,  L=2I+Lk  is  the  total  length). 

Let  us  apply  the  linearized  equation  of  the  flow  in  the  cylindrical  coordinate  system  [5] 


(l  - M 2 )■ 


32cj> 

dx 2 


32<J>  1 3<f> 

dr 2 r dr 


(2.1) 


where  = <I>  * / V0Rn  is  the  dimensionless  flow  velocity  potential, 

<1>*  is  the  flow  velocity  potential, 

Vo  is  the  free  stream  velocity, 

M=  Vo/ao  is  the  Mach  number, 
do  is  the  free  stream  speed  of  sound 

We  assume  that  the  parameter  of  the  cavity  thinness  has  the  same  order  as  the  parameter  of  the  cone- 
cavitator  thinness.  The  parameter  of  thinness  e equals  the  ratio  of  the  cone  radius  to  its  altitude 
(e=l/l=tga).  Let  us  represent  the  potential  O in  the  following  way:  &=(p+x,  where  (p  is  the 
dimensionless  perturbed  velocity  potential  which  is  aimed  at  zero  when  x and  r are  aimed  at  infinity. 
The  dynamical  boundary  condition  on  the  cavity  surface  is  written  as  follows  [6,7]  (we  neglect  the 
members  that  have  the  order  of  trifle  exceeding  £4ln~£ ) 


2 


dip 

dx 


dr 

v / 


= <7 


(2.2) 
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For  the  slender  axisymmetrical  bodies  the  perturbed  velocity  potential  is  defined  by  the  method  of 
sources  and  sinks  distributed  along  the  axis  of  symmetry.  For  the  subsonic  flow  the  potential  (p  satisfying 
the  equation  (2.1)  has  the  form  [5] 


Aj 

r q{x  i )dx  i 

4 n j 

i 

1 V 

t(x  - X,  ) + 1 

(l  - M 2 ) 

1 r2 

(2.3) 


where  q (x j)  is  the  intensity  of  sources  and  sinks  on  the  axis  of  symmetry.  Near  the  axis  of  symmetry  the 
perturbed  velocity  potential  is  represented  by  the  following  asymptotic  equation  [6,7] 


(p  = ^-S'{x)\nr  + g(x) 
in 


(2.4) 


where  S'(x)  is  the  derivative  of  the  dimensionless  area  of  the  slender  body  cross  section  with  respect  to  x 
coordinate. 

The  logarithmic  potential  is  the  main  part  of  expression  (2.4),  and  some  function  g(x)  is  added  to  it.  The 
function  g(x)  is  determined  from  the  condition  of  matching  with  the  potential  of  the  distributed  sources 
(2.3)  since  the  logarithmic  potential  does  not  satisfy  the  condition  at  infinity.  Also,  the  function  g(x) 
takes  into  account  the  influence  of  compressibility.  After  substitution  g(x)  we  obtain  the  following 
expression  for  the  perturbed  velocity  potential  near  the  surface  of  the  body  and  the  cavity  [7] 


(p  = — S'(x)ln 
4 n 


(l  - M2)r2 
4 x(L  - x ) 


1 'pS  'tv,  )-£'(*) 

4 *o  |*“*.| 


(2.5) 


Substituting  equation  (2.5)  in  (2.2)  we  obtain  the  integro-differential  equation  for  the  cavity  profile  [8] 


4—  + „ a !» I1  zgL  ih  _ f 

2 u 4 x(L  - x)  •* 


u r n”(xx  )-  u"(x) 


tlx , — 


-7  <,7l)-‘|('t)Ai - J - 2CT 

J r — y J 


; \x~x\\ 


,+Lk 


(2.6) 


U = R2  , M , = 


Rp  ,u'(0)=0,u'(L)=0 
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The  cavity  profile  is  determined  by  altering  its  radius  along  the  coordinate  x:  R(x).  Along  the  coordinate 
x the  radius  R]  of  the  cone-cavitator  and  the  closing  cone  respectively  change  to  R/=ex  and  R/=£(L-x). 
The  following  boundary  conditions  are  added  to  the  integro-differential  equation  (2.6) 


x = I:R  = l,R'  = e 
x — / + L t . R = \ , R = —£ 


We  seek  the  solution  for  the  whole  area  of  the  cavity  by  expanding  the  asymptotic  rows  with  the  small 
parameter  £. 


tadr'l 

+R2J\n\)  +... 

6 J 

{ £ ) J 

fi  0 

In— 

+(j0+<j{\n\ 

l £ J 

{ £ ) 

(2.7) 


After  substitution  of  the  rows  (2.7)  into  (2.6)  and  conservation  of  the  two  members  of  the  rows  the 
integro-differential  equation  (2.6)  is  transfonned  to  two  differential  equations.  The  first  equation  is 
obtained  from  the  equality  of  the  members  at  £4lnl/£ , the  second  one  follows  from  the  equality  of  the 
members  at  £4 . The  first  differential  equation  with  the  boundary  conditions  is  written  as  follows 


^ = -2a 

dx2  1 


R2{l)=l2,R2(l  + Lk)  = l2 


dR2 

dx 


X-l 


= 21 


(2.8) 


The  solution  of  the  equation  (2.8)  has  the  form 


R l =a](a-x\x-b) 


21 


L Zr  IL 


=—’a=-+\\- — — - 


4 2 


L IL 


4 2 


(2.9) 


After  reduction  of  equations  (2.7)  and  (2.9)  to  the  dimensional  form  and  some  transformations  the  first 
approximation  of  the  cavity  shape  is  written  as  follows 


R 


2 R 


2 / 


IL. 


2 , IL 

■ x + xL 


) 


(2.10) 
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If  we  transfer  the  origin  of  the  coordinate  system  from  the  apex  of  the  cone  to  its  base  then  the  cavity 
profile  takes  the  form  of 


R 2 = 


2K 

IL. 


+ xL  i + - 


IU 


(2.11) 


For  the  stationary  cavity  we  can  write  the  coordinate  x and  the  length  Lk  as  following: 
x = V0t,Lk  — 2 V0(t,  where  t is  the  time.  Substituting  x and  L * in  (2.11)  and  taking  into  account  that 

when  t=tk  the  equality  R=Rk  is  satisfied,  where  /?/,  is  the  cavity  mid-section  radius,  we  obtain  the 
following  equation 


R2  -R]  t f 

Rl  ~ Rl  ~ h k 


(2.12) 


Equation  (2.12)  corresponds  with  equation  (1.18)  if  S0  = kR]  . Thus,  the  principle  of  independence  of 

the  cavity  sections  expansion  and  the  slender  body  theory  give  the  same  result  - we  can  consider  an 
ellipsoid  of  revolution  to  be  the  first  approximation  of  the  shape  of  the  stationary  cavity.  Taking  into 
account  that  the  first  approximation  is  obtained  for  the  compressible  fluid  we  can  draw  a more  general 
conclusion  - the  principle  of  independence  of  the  cavity  sections  expansion  is  satisfied  in  the  subsonic 
flow  of  the  compressible  fluid.  For  the  case  of  the  non-stationary  cavities  in  the  incompressible  fluid  the 
conclusion  of  the  principle  of  independence  from  the  slender  body  theory  has  been  performed  in  [9] 
where  the  so-called  “circular  model”  was  applied. 

Let  us  consider  the  results  obtained  for  the  stationary  cavities  by  the  numerical  methods.  We  can  divide 
the  numerical  methods  into  three  groups.  The  first  group  can  be  called  ‘the  boundary  integral  equation 
method’.  In  this  method  the  stream  function  is  applied  and  the  surfaces  of  a cavity  and  a cavitator  are 
represented  by  the  vortex  layer.  The  cavitation  flow  is  found  with  the  help  of  the  numerical  solution  of 
the  vortex  layer  integral  equations  [10,11,12].  The  boundary  element  method  based  on  the  Green’s 
integral  is  related  to  the  second  group.  In  this  method  on  the  cavitator  surface  the  sought  function  is  the 
perturbed  velocity  potential;  on  the  cavity  surface  the  sought  function  is  the  normal  derivative  of  the 
perturbed  velocity  potential  [13,14].  The  finite-difference  method  represents  the  third  group  [15,16]. 

The  published  works  on  the  numerical  calculations  of  the  stationary  cavitation  flows  contain  the  results 
on  the  cavitation  drag  coefficient  of  a cavitator,  the  cavity  mid-section  radius  and  the  cavity  length.  In 
these  works  the  difference  in  the  results  of  the  numerical  calculations  does  not  exceed  several  percent. 
The  numerical  calculations  have  shown  that  the  stationary  cavity  shape  is  close  to  an  ellipsoid  of 
revolution  for  the  cavitation  numbers  that  have  the  order  lO’^lO"1.  This  result  agrees  with  the  principle 
of  independence  of  the  cavity  sections  expansion. 

As  an  example,  let  us  consider  the  results  of  the  numerical  calculations  obtained  by  the  author  [16]. 
Using  the  finite-difference  method  the  author  has  calculated  the  cavities  past  a disk  in  the  subsonic  flow 
of  the  compressible  fluid.  The  Riabouchinsky  scheme  was  applied  and  the  Mach  numbers  were  located 
on  the  range  0<M<0.95.  The  numerical  calculation  has  shown  that  on  the  whole  range  of  the  Mach 
numbers  the  cavity  profile  is  close  to  an  ellipsoid  of  revolution  which  is  represented  by  the  equation 
(1.18).  An  ellipsoid  joins  with  the  cavity  profile  in  the  vicinity  of  the  disk;  as  the  point  of  contact  we  can 
take  the  point  with  the  coordinates:  x=2R„,  R=2  Rn.  Figure  4 illustrates  the  cavity  profile  in  the  vicinity 
of  the  disk  (the  geometric  dimensions  are  scaled  by  the  radius  of  the  disk  Rn). 
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Fig.  4 The  front  of  the  caivty 


In  Fig.  4 the  continuous  curve  represents  the  result  of  the  numerical  calculation  by  the  boundary  integral 
equation  method  for  the  incompressible  fluid  [10],  the  author’s  numerical  results  in  the  subsonic  flow  for 
the  different  Mach  numbers  are  shown  by  the  points  [16].  Figure  5 shows  the  cavity  profiles  in  the 
compressible  fluid  (M=0.8)  and  in  the  incompressible  fluid  (M=  0)  for  the  same  cavitation  number 
(7=0.0235  [16]. 


0-OM235  M=Q.  8 


Both  of  the  cavity  profiles  are  close  to  an  ellipsoid  of  revolution.  The  dimensions  of  an  ellipsoid  for 
M=  0.8  exceed  the  analogous  dimensions  for  M=  0 since  the  cavitation  drag  coefficient  for  M=  0.8  exceeds 
that  for  M=0. 
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3.  The  application  of  the  equation  of  conservation  of  energy  for  determining  the 
stationary  cavity  shape  at  cr— >0 


The  equation  of  conservation  of  energy  on  the  given  cavity  section  (1.1)  is  equivalent  to  the  principle  of 
independence  of  the  cavity  sections  expansion  since  all  the  three  quantities  that  construct  equation  (1.1) 
are  independent  from  both  the  following  and  the  previous  cavitator  motion.  The  sum  of  the  kinetic  and 
potential  energies  on  the  section  X (Fig.  1)  is  determined  by  the  value  of  the  cavitator  drag  W at  the 
moment  when  a cavitator  passes  the  section  X.  The  law  of  the  cavity  section  expansion  (1.15)  obtained 
from  the  energy  equation  ( 1 . 1 ) is  independent  from  the  following  and  the  previous  cavitator  motion.  This 
law  represents  the  principle  of  independence.  It  should  be  noted  that  the  energy  equation  similar  to 
equation  (1. 1)  was  applied  in  [17]  in  order  to  estimate  the  axisymmetric  cavity  shape.  However,  in  [17] 
the  kinetic  energy  in  the  wake  was  determined  for  the  radial  flow  on  the  plane  circular  layer.  In  chapter  1 
in  contrast  to  [17]  the  kinetic  energy  (formula  (1.2))  is  determined  along  the  stream  tube  leaning  on  the 
expanding  cavity  hole. 

The  principle  of  independence  is  some  approximation  to  the  reality.  However,  a large  number  of 
experiments  have  confirmed  its  accuracy  for  both  the  stationary  and  non-stationary  cavities.  The 
principle  of  independence  agrees  with  the  results  obtained  from  the  slender  body  theory  for  the  stationary 
and  non-stationary  cavitation  flows.  Furthermore,  the  calculation  results  of  the  stationary  cavities  using 
the  numerical  methods  show  that  the  cavity  shape  is  close  to  an  ellipsoid  of  revolution  for  the  cavitation 
number  that  have  the  order  10'2-s-10"'.  Also,  the  ellipsoidal  form  of  the  stationary  cavity  has  been 
obtained  from  the  principle  of  independence  as  we  neglect  the  value  of  v‘/2  compared  to  the  value  of 
AP/p  in  equation  (1.5).  Let  us  prove  that  the  principle  of  independence  of  the  cavity  sections  expansion 
or  the  energy  equation  (1.1)  are  applicable  for  the  special  case.  That  is  the  cavitation  number  is  equal  to 
zero  (the  infinite  cavity  length). 

When  the  cavitation  number  equals  to  zero  then  AP  = - Pk  = 0 and  the  potential  energy  in  the  wake 

E (equation  (1.1))  also  equals  zero.  The  kinetic  energy  in  the  wake  is  determined  by  formula  (1.2).  Let 
us  determine  the  quantities  of  the  velocities  on  the  cavity  surface  at  the  zero  cavitation  number  (Fig.  6). 


y 


x 


Fig.  6 The  absolute  velocity  of  the  fluid  particle  on  the  cavity  surface  at  <7=0 


The  absolute  velocity  of  the  fluid  particles  on  the  cavity  surface  v is  determined  as  the  vector  sum  of  the 
transport  velocity  Vo  directed  along  the  x axis  and  the  relative  velocity  directed  along  the  tangent  to  the 
cavity  surface  (for  the  zero  cavitation  number  the  relative  velocity  is  equal  to  Vo).  In  Fig.  6 the  section  of 
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the  cavity  surface  is  designated  as  s-s.  From  Fig.  6 we  obtain  the  following  correlations  for  the  velocity 
components. 


_r,  • OC  ' . 3(0  ,r  . 

v = 2L,sm— , v ,=  R = Vn since,  v = — = Ln  since 

0 2 y 0 " dn  (3.1) 

vn=vy=R,  v~R 


where  ce  is  the  angle  of  inclination  of  the  cavity  surface  to  x axis.  Correlations  (3.1)  are  approximately 
satisfied  for  the  small  cavitation  numbers  different  from  zero.  They  were  used  in  chapter  1.  For  the  zero 
cavitation  number  the  energy  equation  (1.1)  we  can  write  in  the  following  form  (taking  into  account 
equation  (1.2)  and  correlations  (3.1)) 

- p(pnRR  = W (3.2) 


Let  us  substitute  (p  from  equation  (3.2)  in  the  boundary  condition  (1.5).  We  assume  that  the  approximate 
equality  v ~ R is  satisfied  as  O' is  a small  quantity  (correlations  (3.1)).  As  a result  we  obtain  the  equation 


W_d_ 
Kp  dt 


RR 


\ — = 0 
2 


(3.3) 


The  first  integral  of  the  differential  equation  (3.3)  has  the  form  [1] 


R=1  W 

X‘\np(lnR/R„  + A) 


(3.4) 


where  A is  some  constant. 

For  the  stationary  cavity  the  equality  dx=Vdt  is  satisfied.  Let  us  introduce  the  following  designations: 
X*=x/Rn,  R*=R/Rn.  Equation  (3.4)  is  written  as  follows 

dR* 
dx* 


2 R*J\nR*+A 


(3.5) 


Let  us  join  the  cavity  profile  (3.5)  with  the  cavity  profile  in  the  vicinity  of  the  disk  (1.19).  As  the  point  of 
contact  we  take  the  point  with  the  coordinates  x*=2,  R*=  1.92.  At  the  point  of  contact  we  determine  the 
derivative  dRVdx * from  equation  (1.19)  and  substitute  it  in  equation  (3.5).  As  a result  we  determine  the 
constant  A.  According  to  the  data  presented  in  [1]  the  constant  A equals  0.845.  We  obtain  the  cavity 
profile  from  equation  (3.5)  in  the  form  of 


-v/lnw  + Adu 


x*  = 2 + 


(3.6) 
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where  u is  the  variable  of  integration.  The  asymptotic  law  of  the  axisymmetric  cavity  expansion  has  been 
obtained  for  the  zero  cavitation  number  in  [18.19],  This  law  can  be  presented  as  follows 


R*~4j4C~x 


(3.7) 


In  [20]  the  law  of  the  cavity  expansion  was  generalized  for  the  case  of  subsonic  compressible  fluid  flow. 
In  [20]  it  is  shown  that  equation  (3.7)  has  the  same  form  for  both  incompressible  and  compressible  fluid. 
The  compressibility  exerts  influence  on  the  law  of  the  cavity  expansion  by  means  of  the  drag  coefficient 
Cx,  which  is  dependent  upon  the  Mach  number. 

Let  us  compare  the  cavity  profiles  past  a disk  in  incompressible  fluid  for  the  zero  cavitation  number 
(Cx=0.82).  We  compare  the  profile  obtained  from  the  law  of  conservation  of  energy  in  the  wake 
(equation  (3.6))  with  the  profile  expressed  by  the  asymptotic  law  (3.7).  The  results  of  this  comparison 
are  shown  in  the  table. 


Table 


X* 

R*  eq.  (3.7) 

R*  eq.  (3.6) 

5 

2.67 

2.59 

10 

3.45 

3.36 

25 

5.02 

4.90 

50 

6.77 

6.63 

10 

9.19 

9.03 

0 

25 

13.88 

13.70 

0 

50 

19.06 

18.86 

0 

10 

26.25 

26.04 

00 

15 

31.69 

31.48 

00 
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It  is  evident  from  the  table  that  the  difference  in  the  cavity  profiles  does  not  exceed  several  percent.  We 
can  conclude  here  that  for  the  zero  cavitation  number  the  results  obtained  based  on  the  energy  equation  in 
the  wake  (the  principle  of  independence)  agree  with  the  well-known  Levinson-Gurevich  asymptotic  (3.7). 

In  chapters  1 and  3 we  have  determined  the  cavitation  profiles  based  on  the  energy  equation  (the 
principle  of  independence)  and  have  applied  the  dynamical  boundary  condition  (1.5)  for  two  cases: 

1.  when  we  can  neglect  the  value  of  v2 /2  compared  to  the  value  of  AP/p  in  the  middle  section  of  the 
cavity 

2.  when  the  value  AP equals  zero  (the  zero  cavitation  number). 


In  the  first  case  the  cavity  shape  agrees  with  the  experimental  data  and  the  theoretical  results  for  the 
cavitation  numbers  that  have  the  order  lO^lO"1.  In  the  second  case  the  cavity  shape  agrees  with  the 
Levinson-Gurevich  asymptotic  at  <7=0.  Thus,  we  do  not  investigate  the  cavities  on  the  range  0«7<0.005. 
Let  us  extend  the  application  of  the  principle  of  independence  to  this  range  of  the  cavitation  numbers. 

Recently  G.  V.  Logvinovich  has  considered  the  general  case  - two  members  v'/2  and  AP/p  are  taken  into 
account  in  the  dynamical  boundary  condition  (1.5)  [21],  Let  us  determine  the  stationary  cavity  shape  for 
this  case.  The  potential  energy  is  defined  as  E=7lR2 AP  and  equation  ( 1 .4)  is  written  in  the  following 
form 


- ptpnRR  + kR2 AP  = W 


(3.8) 


After  substitution  of  (p  from  equation  (3.8)  into  the  dynamical  boundary  condition  (1.5)  using  the 
approximate  equality  v ~ R we  obtain  the  following  equation 


W d 

2 1 2 

AP  d 

' R 2 

Ap 

up  dt 

A, 

p dt 

A 

P 

(3.9) 


As  a result  of  the  transfonnations  the  first  integral  of  differential  equation  (3.9)  for  the  stationary  cavity 
can  be  written  as  follows  [21] 


dR*  [C7j_  ll-oR*2  !CX 
dx*  V 2 R * V In R * +A 


(3.10) 


At  0=0  differential  equation  (3.10)  is  transfonned  to  equation  (3.5).  Using  equation  (3.10)  the  cavity 
profiles  were  calculated  for  the  wide  range  of  cavitation  numbers  [21].  For  the  cavitation  numbers  that 
have  the  order  10'2  the  calculations  on  the  basis  of  equation  (3.10)  have  shown  that  the  cavity  shape  is 
close  to  an  ellipsoid  of  revolution.  This  result  agrees  with  the  well-known  experimental  and  theoretical 
results. 

Thus,  on  the  basis  of  the  carried  out  investigations  we  can  conclude  that  the  principle  of  independence  of 
the  cavity  sections  expansion  deduced  from  the  equation  of  conservation  of  energy  in  the  wake  is 
applicable  to  the  wide  range  of  the  cavitation  numbers  0<O<0.1.  On  this  range  the  stationary  cavity 
profiles  can  be  calculated  by  the  single  formula  (3.10). 
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4.  The  calculations  of  the  cavities  in  the  gravitational  field  on  the  basis  of  the 
principle  of  independence  of  the  cavity  sections  expansion 


The  principle  of  independence  of  the  cavity  sections  expansion  makes  the  calculation  of  the  cavity  profile 
easier  in  the  arbitrary  variable  pressure  fields.  Based  on  this  principle  the  cavity  shape  is  determined  if 
we  know  the  law  of  the  cavitator  motion  and  the  values  of  the  pressure  at  infinity  and  within  the  cavity. 
Let  us  apply  the  principle  of  independence  for  the  calculation  of  the  vertical  cavities  in  the  gravitational 
field  [22], 


Let  us  consider  the  vertical  axisymmetrical  cavity  formed  by  the  stream  around  the  motionless  disk- 
cavitator  in  the  gravitational  field.  The  stream  can  be  descending  or  ascending,  the  value  of  the  stream 
velocity  V is  the  constant  (Fig.  7). 


Fig.  7 The  coordinate  systems  for  the  calculation  of  the  cavity  profile  in  the  gravitational  field 


At  the  level  of  the  cavitator  the  free  stream  pressure  P„  (the  pressure  at  infinity)  is  the  constant.  Let  us 
introduce  the  vertical  coordinate  axis  x directed  along  the  stream  velocity  and  the  symmetry  axis  of  the 
cavity  from  the  cavitator  center.  On  the  section  with  coordinate  x the  free  stream  pressure  (the  pressure 
at  infinity)  is  expressed  as  following 

PSx)=Pn  ± Pg-X  (4-1) 


where  g is  the  gravitational  acceleration,  the  sign  ‘plus’  corresponds  with  the  descending  stream,  the  sign 
‘minus’  corresponds  to  the  ascending  stream. 

We  determine  the  profile  of  the  stationary  vertical  cavity.  For  the  stationary  cavity  the  pressure  within  it 
is  the  constant  - PK— const.  The  area  of  the  cavity  cross  section  fixed  relatively  the  motionless  fluid 
changes  according  to  equation  (1.15).  Integral  of  equation  (1.15)  is  equation  (1.17)  which  describes  the 
change  of  the  area  of  the  cavity  fixed  cross  section.  In  this  example  we  consider  the  cavitator  to  be 
motionless  and  observe  the  cavity  section  which  has  been  formed  at  the  edge  of  the  cavitator  at  the 
moment  t= 0 and  moves  at  the  velocity  V along  with  the  stream.  The  area  change  of  this  section  can  be 
represented  by  the  following  expression 

it  u 

— Ji/m  | AP(v)c(v 

P 0 0 


S(t)=S0+S0t 


(4.2) 
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where  v and  u are  the  variables  of  integration.  So  is  the  section  initial  area  equal  to  the  cavitator  area.  Let 
us  substitute  the  expression  of  the  external  pressure  (4.1)  in  dependence  (4.2).  Note  that  in  (4.1)  x is 
written  as  x(v)=Vv.  After  integration  of  equation  (4.2)  we  obtain 


kt2  ( 1 

S(t)=S0+S0t-—  AP0  ± — pgVt 

2 P { 3 


Expression  (4.3)  contains  the  difference  APo  of  the  pressure  at  infinity  and  the  pressure  within  the  cavity 
AP0  = Pn  - Vk  . This  difference  is  related  to  the  section  at  the  level  of  the  cavitator. 

Let  us  consider  the  case  of  the  descending  stream.  For  the  stationary  cavity  the  variable  x is  written  as 
x—Vt.  After  substitution  of  (1.14)  into  (4.3)  and  certain  transformations  we  obtain  the  expression  for  the 
vertical  cavity  profile  in  the  dimensionless  form 


2 C C r*2 

r*2  + 


where  R*,  x*  are  the  dimensionless  cavity  radius  and  coordinate  x scaled  by  the  disk  radius  Rn, 
O0  = 2AP0  / pV 2 is  the  cavitation  number  at  the  level  of  the  cavitator, 

Fr  = V / j2Rng  is  the  Froude  number. 


Fig.  8 The  vertical  cavity  profiles  in  the  descending  stream 
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As  an  example,  in  Fig.  8 the  cavity  profiles  calculated  by  formula  (4.4)  are  represented  by  the  continuous 
curves  for  the  three  values  of  the  cavitation  number  a0:  -0.1;  0;  0.1.  The  Froude  number  equals  3.2  and 
it  remains  the  same  for  all  the  three  cavities,  the  constant  a equals  1.5.  The  Levinson-Gurevich 
asymptotic  (3.7)  for  the  zero  cavitation  number  and  weightless  fluid  ( Fr=°° ) is  represented  by  the  broken 
curve. 

Let  us  determine  the  length  of  the  stationary  vertical  cavity  for  the  Riabouchinsky  scheme  (Fig.  7). 
We  assume  that  the  cavity  is  closed  by  the  disk  of  the  same  dimensions  as  the  disk  cavitator,  i.e.  i?*=l  at 
x*=Lii/R„.  From  expression  (4.4)  we  obtain  the  quadratic  equation  for  determining  the  cavity  length 


( 


V 


\2 

) 


3 ciFr2 


■ + 


hi 


2 = 0 


(4.5) 


The  solution  of  equation  (4.5)  can  be  represented  in  the  following  form 


2 R,Fr 


-auFr  + 


o:Fr 2 + — 


3 a 
2 


(4.6) 


In  Fig.  9 dependence  (4.6)  is  represented  by  the  continuous  curve,  the  points  show  the  experimental  data 
obtained  for  different  values  of  the  Froude  and  the  cavitation  numbers  [22].  Figure  9 shows  good 
agreement  between  the  theoretical  dependence  (4.6)  and  the  experimental  data. 


Fig.  9 The  comparison  of  the  calculation  of  the  vertical  cavity  length  based  on  formula  (4.6) 
with  the  experimental  data 


Let  us  apply  the  principle  of  independence  for  the  investigation  of  the  cavity  deep  closure  which  occurs 
at  the  water  entry  of  bodies  [23].  At  the  moment  of  intersection  of  the  water  surface  the  cavity  is  formed 
past  a body.  This  cavity  communicates  with  the  atmosphere,  consequently  in  the  initial  period  of 
submersion  the  gas  pressure  within  the  cavity  is  close  to  atmospheric  pressure  and  the  cavitation  number 
is  close  to  zero.  However,  as  the  body  sinks  the  external  hydrostatic  pressure  increases  and  under  the 
influence  of  the  excess  pressure  the  cavity  section  closes.  The  cavity  detaches  from  the  atmosphere  and 
the  so-called  deep  closure  occurs  (Fig.  10). 
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Fig.  10  The  cavity  deep  closure  at  the  water  entry  of  a cavitator 


Let  us  consider  the  inclined  rectilinear  water  entry  of  the  cavitator  at  the  constant  velocity  V.  We 
designate  the  entry  angle  as  0 (the  angle  of  inclination  of  the  trajectory  to  the  free  surface).  The  time  t is 
measured  from  the  instant  of  the  cavitator’ s initial  contact  with  the  water  surface.  We  observe  the  section 
which  has  been  formed  at  the  moment  t\  on  the  depth  H\  (Fig.  1 1). 


V 


Fig.  11  The  scheme  for  the  calculation  of  the  deep  closure 


We  consider  the  pressure  within  the  cavity  equal  to  the  atmospheric  pressure  then  the  excess  pressure  on 
the  considered  section  is  the  constant  and  can  be  written  in  the  form  of 


AP  = pgHx  = pgL,  sin  6 


(4.7) 
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The  equation  of  the  cavity  section  expansion  is  described  by  equation  (1.15).  The  integral  of  equation 
( 1 . 1 5)  is  written  as  follows 


S(Hl  ,t)=S0+S0{t-tl)-—-{t-tl)2  (4.8) 

2 P 

where  So  is  the  cross  section  initial  area  equal  to  the  cavitator  area.  From  equations  (4.7)  and  (4.8)  we 
can  determine  the  cavity  profile  at  the  water  entry  of  the  disk-cavitator  if  we  apply  the  obvious 
correlations  - x=V(t  — /,);  L{=L  — x where  x is  the  coordinate  measured  from  the  cavitator  center 
along  the  trajectory  of  the  cavitator,  L is  the  cavitator  path  passed  after  the  intersection  of  the  water 
surface. 


S(x)=S0  + y-x' 


kgL  sin  0 2 ^ kg  sin  0 3 


2V1 


2V 2 


(4.9) 


We  determine  the  cavity  length  from  the  condition  that  the  cavity  is  closed  by  the  disk  with  the  same 
dimensions  as  the  disk-cavitator.  From  (4.9)  we  obtain  the  quadratic  equation  for  determining  the  cavity 
length 


kg  sin  0 
2V2 


kgL  sin  0 T S0 

2V2  V 


(4.10) 


The  solution  of  equation  (4. 10)  is  written  in  the  form  of 

L =---]  1 

2 2 V kgL2  sin  0 


(4.11) 


It  is  obvious  that  the  maximum  cavity  length  corresponds  to  the  deep  closure.  It  follows  from  (4.1 1)  that 
the  value  of  Z./{  is  maximum  when  the  radicand  expression  equals  zero.  Then  the  value  of  Li<  equals  LI 2, 
i.e.  the  cavity  closes  halfway  passed  by  the  cavitator.  The  moment  of  the  deep  closure  is  determined 
from  the  equation 


8S0V 

kgL2  sin0 


(4.12) 


After  substitution  of  correlations  (1.14)  and  the  equality  L=Vtd  in  (4.12),  for  the  dimensionless  deep 
closure  time  we  obtain  the  expression 


2a 

sin0 


C = 2 Fr. 


(4.13) 
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where  td  — Vtd  / Rn  is  the  dimensionless  deep  closure  time,  is  the  deep  closure  time. 

Previously  expression  (4.13)  was  obtained  in  a different  way  in  [24].  The  path  passed  by  the  cavitator  at 
the  moment  of  the  deep  closure  is  determined  from  the  equality  L = Rjd . 


5.  The  application  of  the  principle  of  independence  for  the  calculation  of  the  non- 
stationary cavities  with  variable  internal  pressure 


The  area  of  the  cavity  cross  section  fixed  relatively  to  the  motionless  fluid  changes  according  to  equation 
(1.15).  Equation  (1.15)  expresses  the  principle  of  independence.  The  integral  of  equation  (1.15)  is 
equation  (1.17)  that  describes  the  area  change  of  the  cavity  fixed  cross  section.  Let  us  go  on  from 
equation  (1.17)  to  the  expression  that  represents  the  cavity  profile  at  the  fixed  time  [25],  In  order  to  do 
that  in  equation  ( 1 . 1 7)  we  introduce  the  time  t which  is  common  for  all  the  sections.  At  the  moment  t the 
cavitator  has  the  coordinate  H.  Let  us  designate  the  moment  of  the  section  formation  with  the  coordinate 
h as  t\  (t\=t\(h)<t=t\(H)).  Then  equation  (1.17)  transforms  to 


j i u 

S(h,t)=  S0+S0(t-tl)-—jdujAP(h,v)d 


Pi, 


(5.1) 


where  So  is  the  cavitator  area,  S0  is  the  constant  that  defines  the  initial  velocity  of  the  cavity  expansion, 

v and  u are  the  variables  of  integration.  At  the  fixed  moment  t equation  (5. 1)  expresses  the  cavity  profile 
S(h).  At  the  fixed  moment  t\,  i.e.  /t=const,  equation  (5.1)  describes  the  expansion  of  the  fixed  cavity 
cross  section.  We  can  transform  the  double  integral  in  (5. 1 ) to  the  line  one  as  follows.  Let  us  change  the 
order  of  integration  in  the  double  integral 

t U t t t 

J du  J AP(h,v)dv  = J AP  (/t,  v)jv  j du  = J (/  - v)\P  (h,  v)dv  (5.2) 

/|  f|  /j  v 

After  substitution  of  (5.2)  in  (5. 1)  we  obtain 

7 t 

S(h,t)=  S 0 + S0  (t  - t,  ) [ (/  - v)VP  (h,v)dv  (5.3) 

pi 


Equation  (5.3)  is  especially  convenient  for  the  numerical  calculation  of  the  cavity  profile  when  the 
pressure  within  it  changes  arbitrarily.  As  an  example  we  consider  the  calculation  of  the  cavity  profile 
when  the  gas  pressure  within  it  changes  according  to  the  harmonic  oscillations  [26].  The  time 
dependence  of  the  gas  pressure  is  written  as  follows 


Pk(t)=Pk  + A sin  cot 


(5.4) 


8-22 


where  Pk  is  some  constant  pressure,  CO  is  the  circular  frequency,  A is  the  amplitude  of  oscillations.  We 

consider  the  external  pressure  P,„(h,t)  to  be  a constant  and  equal  to  Poo.  After  substitution  of  (5.4)  in 
(5.3)  or  (5.1)  we  obtain  the  profile  of  the  pulsating  cavity 


S(x,t)=S0  +-j^x~ 


MP, 


2 pV2 


0 x2 


kA 

par 


\ 

sin[ft>  t + <p(x)] 


(5.5) 


where  x = V{t-t])  is  the  coordinate  measured  from  the  cavitator  along  the  stream  velocity, 
A P0  = P„  — Pk  is  the  stationary  excess  pressure 


F2 


cox 

-cos  — 
V 


cox  _ . cox 

2 sin  — 

V V 


(5.6) 


The  first  three  members  in  (5.5)  express  the  profile  of  the  stationary  cavity,  the  last  member  describes  the 
perturbation  of  the  cavity  profile.  The  amplitude  of  the  perturbation  of  the  cavity  profile  increases  to  the 
downstream  end  of  the  cavity.  The  increase  of  the  perturbations  to  the  downstream  end  of  the  cavity  is 
illustrated  by  Fig.  12  which  shows  the  dependence  (5.6). 


Fig.  12  The  increase  of  the  perturbations  to  the  downstream  end  of  the  cavity 


Also,  the  increase  of  the  perturbations  to  the  downstream  end  of  the  cavity  is  confirmed  by  the 
experiments  performed  with  the  pulsating  cavities  [26]. 

If  beforehand  the  gas  pressure  within  the  cavity  is  not  known  then  it  is  necessary  to  add  the  subsidiary 
equations  to  equation  (5.3)  or  (1.15).  The  equation  for  calculating  the  gas  volume  within  the  cavity  is 
one  of  them 


<2(0  = 


*«(») 

Jb 


H (/) 


(h , t )-  S h (h , / )\dh 


(5.7) 


where  Q(t)  is  the  volume  contained  between  the  cavity  boundary  and  the  body  within  the  cavity, 
St  (h,  t)  is  the  area  of  the  body  cross  section; 

H(i)  is  the  coordinate  of  the  cavitator, 
hc  (t)  is  the  coordinate  of  the  cavity  closure. 


8-23 


The  equation  of  the  gas  mass  change  within  the  cavity  has  the  form  of 


dm(t) 

— = m,  - m, 

dt  s ‘ 


(5.8) 


m(t)  is  the  gas  mass  within  the  cavity; 
ths  is  the  mass  gas-supply  rate  to  the  cavity; 

ml  is  the  mass-leakage  rate  from  the  cavity. 

The  following  is  the  equation  of  the  gas  state  within  the  cavity 


pk(t)  Q(t)  = RT(t) 
m{t) 


where  R is  the  gas  constant, 

T(t)  is  the  gas  temperature  within  the  cavity.  This  temperature  is  dependent  upon  the  thermodynamic 
process.  In  most  cases  we  can  consider  the  thermodynamic  processes  within  the  cavity  to  be  isothermal 
and  assume  that  the  gas  temperature  within  the  cavity  is  equal  to  the  temperature  of  the  surrounding  fluid. 

In  general  the  equation  of  the  cavitator  motion  appears  like 


dV 

dt 


F{V,H,Pk,hc„. ) 


(5.10) 


where  V is  the  velocity  vector. 

The  system  of  the  equations  (5.3),  (5.7)-(5.10)  describes  the  dynamics  of  the  non-stationary  cavity  in  the 
closed  form.  In  fact,  we  can  solve  this  system  by  the  numerical  method  at  each  time  step  substituting  the 
cavity  profile  for  the  discrete  set  of  the  cross  sections.  Also,  in  this  case  using  the  numerical  method  we 
can  solve  equation  (1.15)  instead  of  equation  (5.3). 

E.  V.  Parishev  [25]  has  transformed  the  integral  correlations  (5.3)  and  (5.7)  in  the  nonlinear  differential 
equations  with  the  lagging  argument.  It  turned  out  that  in  general  the  dynamics  of  the  non-stationary 
cavity  is  described  by  the  system  of  the  nonlinear  differential  equations  of  the  sixth  order  with  the 
variable  lag  (excluding  the  equations  of  the  cavitator  motion).  In  particular  cases  these  equations  are 
simplified.  For  example,  in  case  of  the  axisymmetric  cavity  formed  past  a cavitator  which  is  moving  at 
the  constant  velocity  in  the  weightless  fluid  and  with  the  assumption  that  the  gas  mass  within  the  cavity  is 
the  constant  (gas-supply  and  gas-leakage  are  absent),  the  equation  for  the  small  pressure  oscillations  of 
the  gas  within  the  cavity  has  the  form  [27] 


:Pk(t*)+' 


-Pk{< 


' — T *)=  0 


(5.11) 
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where  Pk  is  the  gas  pressure  within  the  cavity  (small  oscillations  relatively  to  the  equilibrium),  t*=t/T-  is 
the  dimensionless  time,  x*=x/T  is  the  dimensionless  lag,  X=Lk/V  is  the  dimensional  lag,  Lk  is  the  cavity 
length  at  the  equilibrium,  V is  the  cavitator  velocity.  The  time  scale  T is  determined  by  the  expression 


nQoP 

]P*.kVT 


where  Qu  is  the  cavity  volume  at  the  equilibrium,  Pt  is  the  gas  pressure  within  the  cavity  at  the 

equilibrium,  1 /n  is  the  polytropic  index  of  the  thermodynamic  process  within  the  cavity,  for  the 
isothermal  process  n equals  1. 

Equation  (5.11)  has  the  single  basic  parameter  equal  to  the  dimensionless  lag  T*.  As  a result  of 
transforming  equation  (5.12)  this  parameter  can  be  represented  as  following 


T*  = . 


inf Eu_ _ N 


where  Eu  = 2 Px  / pV2  is  the  Euler  number,  <r  is  the  average  cavitation  number  defined  by  the  pressure 

n(1- 

Thus,  all  the  cavity  properties  described  by  equation  (5.1 1)  are  defined  by  the  value  of  the  ratio  Eu/a. 
Also  the  experimental  works  [28,29]  show  that  the  phenomenon  of  the  artificially  inflated  cavity 
pulsation  is  defined  by  the  ratio  a/Eu.  Equation  (5.11)  has  both  the  steady  and  unsteady  (pulsating) 
solutions.  For  the  vapor  cavities  (o~Eii)  and  the  ones  close  to  them  the  solutions  are  steady.  Flowever, 
when  X*  > k42  the  solutions  are  unsteady.  This  result  corresponds  to  the  well-known  features  of  the 
vapor  cavity  which  does  not  pulsate.  Also,  the  experiments  in  [28,29]  show  that  the  cavities  with  a weak 
gas-supply,  i.e.  those  that  are  close  to  the  vapor  ones,  are  steady. 


Fig.  13  The  dependence  of  the  dimensionless  circular  frequency  of  cavity 
pulsation  co * on  the  parameter  T* 


In  Fig.  13  the  dependence  of  the  dimensionless  circular  frequency  of  the  cavity  pulsation  (o*  (a>*=coLi/x* 
V,  where  (o  is  the  dimensional  circular  frequency)  on  x*  is  represented  by  the  discontinuous  serrate 
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function  with  the  average  value  equal  to  1 . Knowing  the  dimensionless  frequency  of  pulsation  we  can 
determine  the  number  of  the  waves  N going  into  the  cavity  length 


CO*  X* 

2k 


where  X is  the  wave  length.  Since  the  dependence  co*(x*)  is  discontinuous  the  dependence  N(x*)  is 
discontinuous  too  and  the  number  of  the  waves  going  into  the  cavity  length  is  close  to  integer.  At  the 
continuous  change  of  r*  the  value  of  N changes  as  the  stepped  function.  The  dependence  N(r*)  is 
represented  in  Fig.  14.  The  points  show  the  experimental  data  [26]. 


61 


Jo  Jo 
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Fig.  14  The  dependence  of  the  number  of  the  waves  going  into  the  cavity  length  N on  the 

parameter  x* 


The  discontinuous  properties  of  the  dependences  co*(T*)  and  N(x*j  agree  with  the  well-known 
experimental  data  [29].  The  experiments  have  shown  [29]  that  the  pulsating  cavities  have  the  regimes 
with  the  fixed  number  of  waves  close  to  integer.  The  experiments  have  also  shown  [29]  that  at  the 
smooth  value  change  of  the  gas-supply  the  number  of  waves  changes  as  the  stepped  function. 

Thus,  the  cavity  instability  is  connected  with  the  gas  elasticity  within  it.  The  underwater  bubble  has  the 
analogous  instability  when  it  is  internally  supplied  by  gas.  We  can  consider  the  cavity  to  be  a dynamic 
system  with  infinite  degrees  of  freedom.  The  cavity  cross  section  changes  influenced  by  the  difference 
between  the  free  stream  pressure  and  the  pressure  within  the  cavity;  the  cavity  volume  changes  as  well. 
The  cavity  volume  change  leads  to  the  pressure  change  within  the  cavity,  etc.  The  cavity  being  a 
dynamic  system  with  these  properties  is  analogous  to  the  vibration  of  a body  on  a spring.  The  coordinate 
of  a body  is  analogous  to  the  cavity  cross  section  area.  The  elastic  force  of  a spring  is  analogous  to  the 
pressure  within  the  cavity.  However,  in  case  of  a cavity  we  have  an  infinite  quantity  of  these  bodies 
(sections)  and  the  pressures  within  the  cavity  is  a function  of  the  whole  cavity  volume,  i.e.  the  cross 
sections  are  connected  by  means  of  the  gas  that  fills  the  cavity.  The  gas-leakage  from  the  cavity  is  the 
factor  of  stability.  When  this  factor  is  active  the  parameter  x*  that  corresponds  to  the  cavity  stability, 
exceeds  the  value  of  njl  . In  addition  to  the  mathematical  model  Parishev  removed  the  gas  along  the 
bound  vortexes  and  obtained  the  accordance  with  Epshtein’s  experimental  data  [26]. 
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Conclusions 


The  principle  of  independence  of  the  cavity  sections  expansion  stated  by  G.  V.  Logvinovich  provides  an 
easier  way  for  investigating  the  non-stationary  cavities  and  the  cavities  with  the  variable  external 
pressure.  This  principle  corresponds  to  the  flow’s  property  around  slender  bodies  and  it  is  stated  as 
following  - each  cavity  cross  section  expands  relatively  to  the  trajectory  of  the  cavitator  center  according 
to  the  certain  law  which  is  dependent  on  the  conditions  at  the  moment  as  a cavitator  passes  the  plane  of 
the  considered  section.  The  cross  section  expands  almost  independently  from  the  following  or  the 
previous  cavitator  motion. 

The  principle  of  independence  is  equivalent  to  the  equation  of  conservation  of  energy  applied  to  the  given 
cavity  section.  The  sum  of  the  kinetic  and  potential  energies  on  the  cavity  section  is  determined  only  by 
the  value  of  the  cavitation  drag  at  the  moment  of  intersection  with  this  section.  The  law  of  the  cavity 
section  expansion  (1.15)  obtained  from  the  energy  equation  (1.1)  is  independent  from  the  following  and 
the  previous  cavitation  motion. 

The  principle  of  independence  is  some  approximation  to  the  reality.  However,  the  numerous  experiments 
have  confirmed  it  being  accurate  for  both  the  stationary  and  non-stationary  cavities.  The  principle  of 
independence  agrees  with  the  results  obtained  from  the  slender  body  theory  for  the  stationary  and  non- 
stationary  cavitation  flows. 

Furthermore,  the  results  of  the  stationary  cavity  calculation  by  the  numerical  methods  demonstrate  that 
the  cavity  shape  is  close  to  an  ellipsoid  of  revolution  for  the  cavitation  numbers  that  have  the  order 
101 2 3 4 5 6 7 8 9-f-10''.  The  ellipsoidal  form  of  the  stationary  cavity  has  also  been  obtained  from  the  principle  of 
independence.  The  principle  of  independence  of  the  cavity  section  expansion  and  the  energy  equation  are 
applicable  for  a special  case  - the  cavitation  number  equals  zero  (the  infinite  cavity  length). 
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